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ABSTRACT 

The statistics of high speed satelhte galaxies, as reported in the recent hterature, can be a powerful 
diagnosis of the depth of the potential well of the host halo, and hence discriminate between competing 
gravitational theories. Naively one expects that high speed satellites are more common in Modified 
Newtonian Dynamics (MOND) than in cold dark matter (CDM) since an isolated MONDian system 
has an infinite potential well, while CDM halos have finite potential wells. In this Letter we report 
on an initial test of this hypothesis in the context of the first generation of cosmological simulations 
utilizing a rigorous MONDian Poisson solver. We find that such high speed encounters are approx- 
imately a factor of four more common in MOND than in the concordance ACDM model of cosmic 
structure formation. 

Subject headings: galaxy: formation - methods: N-body simulations - cosmology: theory - dark 
matter - large scale structure of Universe 



1. INTRODUGTION 



The standard ACDM model fcf..'Koma tsu et al.l[20M l 
explains the formation of cosmological structure in the 
non-linear regime in a hierarchical way, i.e. large struc- 
tures are not formed monolithically b ut by the successiy e 
merging of smaller structures (e.g,. iDavis et"al1 11985D . 
Recent cosmological simulations also support the idea of 
hierar chical formation in MOND gravity ( Llinares ct aL, 
(l200l. but see also analytical models of ISanderd 120081 : 
IZhao et al.ll2008[ ). The hierarchical merging scenario nat- 
urally promotes the picture that we should observe col- 
lisions of galaxies. The question that immediately arises 
is what is the nature of the distribution of the relative 
speed of such encounters. Observationally there is evi- 
dence that some of these collisions actually occur with 
speeds that are not readily re produced by simulations 
of ACDM structure formation d Havashi fc Whitd l2006l : 
iSpringel fc Farrail [2007t iKnebe et al. 2001). There is, 
for example, the famous "Bullet cluster" , an extremely 
high velocity merger between two galaxy clusters whose 
relative speed is between 2500 and 4500 km/sec de- 
pending on the interpretation of the shock speed and 
the method used to infer the collision speed (observa- 
tions/models, a nalytica l /num erical, N-b ody /hydro sim- 
ulations, e.g. iNusseij l2008t ISpringel fc Farrad l2007l : 
lMarkevitchll2006l : IZhadl2007bf ). At first sight, the upper 



limit of this interval is too high and ma y be a problem for 
the standard LambdaGQyi model, but lHavashi fc Whitd 
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(|2006[) showed using the Millenium cosmological simula- 
tion (Springel et al.ll2005i ) that the probability of such an 
event albeit low, is not zero. 

There are a number of such high speed encounters 
in the literature. One example of such a collision is 
the so-called "line-of-sight Bullet", i.e. Abel 576, wit h 
a relative velocity of 3300 km/sec (jPupke et al.ll2007b . 
Furthermore the "Cosmic Train Wreck" Abel 520 is a 
collision with a veloc ity of approximately 1000 km/sec 
(jMahdavi et al.ll2007| ). The "Dark Matter Ring" clus- 
ter C10024-I-1 7 has a specu lated impact velocity of 
3000 km/sec (|Jee et all 12001 and MACS J0025.4-1222 
has two merging components w hose relative velocity was 
measured to be 2000 km/sec (jBradac et al.ll2008( ) . In 
comparison the random dispersion of velocities in these 
clusters is only about 500-1000 km/sec. 

A consequence of any high speed collision of mass con- 
centrations seems to be the decoupling or offsetting of the 
baryons from the dark component. Assuming this being 
the case, additio nal examples o f collisions are given i n 
iJee et all (!2005a'b) but see also lHevmans et al.r (|2008'l . 
This kind of objects, with offsets between baryon and 
DM components, have become what could be considered 
as yet another important standard test that any theory 
for g ravity shou ld pass before being seriously considered 
(e.g. IWilJ ll9931. Simply applying the MOND formula 
to a universe populated only with baryons seems to fail 
this test. Possible solutions could come from many on- 
going efforts to embed the MOND idea in a relativistic 
framework by adding complementary (vector) fields be- 
sides the s tandard Einstein's metric ( BekensteinI 1200 J : 
iZlosnik et al.ii2007i : ^hao 200 7a. 2008) or by the addition 
of neutrinos of various kind lAngus fc McGaughl l2008t 
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iFeix et al.|[200l IZhaoll2008D . However we must have in 
mind, that the same data on the Bullet Cluster would 
have rejected general relativity in its original formula- 
tion without introducing one or more dark matter com- 
ponents plus a cosmological constant. 

The question that arises from all these data is how to 
match th e low probability o f high-speed encounters pre- 
dicted bv lHavashi fc White! ((2006,) for the ACDM model 
with the fact that this type of collisions seems to be com- 
mon in the observable universe. A clue comes from the 
MONDian point of view where the situation seems to be 
more favorable for high velocities. P revious authors have 
noted the deep potential in MO ND (jAngus fc McGaugD 
120081 IZhaoll2007bl : lNusseij|2008l ) is helpful in the context 
of the Bullet cluster. On a smaller scale, high-velocity 
stars have been stu died in the co ntext of the escape speed 
in the Milky Way (jPerets et al.l l2008). It was found that 
MOND can retain stars of higher velocity than CDM, and 
the MONDian escape velocity is more cons istent with the 
RAVE data in the so lar neighbourhood (jFamaev et al.l 
I2007t IWu et al.|[2007[) . Further, the revised (yet still 
discusse d) speed of the Magellanic Clouds also favors 
MOND (jWu et al.l[20M l. The question that previous 
authors cannot address is how to obtain a self-consistent 
strength of the external field in MOND since they lack 
a full cosmological simulation. And as MOND is a non- 
linear theory it violates the strong equivalence principle 
and hence it is mandatory to simulate galaxies within the 
cosmological framework and not in isolation. 

The aim of this paper is not to go further in an explana- 
tion of this kind of systems using MONDian ideas, but to 
study the consequences of a MONDian cosmological toy 
model on the probability of such high speed encounters. 
In order to do this, we study the velocity distribution 
of substructure extracted from cosmological simulations 
that have been run using standard and modified gravity. 
We show that high speed collisions are more frequent in 
MOND than in the concordance ACDM model. 

2. SIMULATIONS 

The analysis presented in this Letter is based upon a 
set of two simulations published in lLlinares et al.l (|2008l ). 
i.e. the ACDM and the OCBMond2 model, respec- 
tively. Both simulations were run in a box with a 
side length of 32/i~^Mpc, using 128"^ particles. They 
were b oth run with a mo dification of the iV-body code 
MLAPM (|Knebe et all 1200 If ). The ACDM model employs 
a background cosmology parametrized by ^dm+b = 0.3, 
f^A = 0.7, and a normalisation of the power spectrum of 
the density perturbation of cts = 0.88. For the MONDian 
simulation, we chose an open universe with neither dark 
matter nor dark energy but characterized by fib = 0.04. 
In order to arrive at a comparable evolutionary stage to 
the ACDM model at redshift z = we had to lower the 
normalisation erg to 0.4 due t o the faster growth of struc- 
tures in MOND (cf.| Saiideri2001; Kne be fc Gibson 2004 ; 



iLlinares et al.ll2008f ). Both simulations were started at 
redshift z = 50 and used a Hubble constant Hq = 70 
km/sec/Mpc. 

We used the MPI version o f the AHF halo fin der^ 
AMIGA's-Halo-Finder (KnoUm ann fc Knebd l2009f) to 
identify objects, which is based on the MHF halo finder of 

^ AHF is freely available from http : //www . aip . de/People/aknebe 



iGill et aT] (|2004j) . For the identification of substructure 
we employed the tool MergerTree that comes with the 
AHF software package. This algorithm was originally de- 
signed to follow halos through time by tracking the mem- 
bership of individual particles, but it can be also used to 
locate the subhalos of a given host. Since particles that 
belong to subhalos will belong also to the corresponding 
host, constructing a merger tree of a halo catalogue with 
itself will provide us with a "subhalo tree" (rather than 
a merger tree). It is important at this moment to make a 
remark about our terminology. We use the term subhalo 
to refer to the largest substructures embedded in host 
haloes. The mass ratio between our host halos and the 
most massive subhalo have a median of 0.23 and 0.15 for 
the MONDian and Newtonian simulations respectively. 
These numbers are in the range of typical mass ratios for 
collisions in mergers of host halos and are well above the 
typical ratio betwe en hosts and real substructures (e.g. 
iMadau et al.ll2008l ). In order to not contaminate our re- 
sult with unvirialized objects we further prune our halo 
catalogue by removing objects with a high virial ratio 
ending up with 64 and 58 objects in the Newtonian and 
MONDian simulations, respectively. 

For more d etails regardin g these simulations, we refer 
the reader to ILlinares et al.l (|2008f ) 

3. ANALYSIS 

While the primary focus of this Letter is the distri- 
bution of the relative velocity of two colliding systems, 
we still need to define a proper normalisation for these 
velocities to correct for the fact that a more massive 
host system will lead to a larger acceleration towards 
its centre. While others referred to the rotational veloc- 
ity at the virial radius of the host for this purpose (e.g., 
iHavashi fc Whit3[200l , we rather use the mass-averaged 
velocity dispersion. 

3.1. Velocity dispersion - Mass relation for MOND and 

CDM 

In the Newtonian case, the velocity dispersion scales 
with the mass M as follows: a oc Vdr = y/GM/R (x 
R oc A/^/"^ cx AfjYn^on' where we used M oc pR^ where 
p is the background density, which depends only on red- 
shift. A similar scaling relation between velocity dis- 
persion a and mass M can be easily obtained for deep 
MOND a oc Vc^r oc {GMbaryonaof/^ OC M.^^o™ for a 
spherical isolated body. Although not rigorous, we find 
that this scaling holds fairly well as a mass-averaged total 
dispersion of the system even in the intermediate MOND 
regime. 

Figure [1] shows the a — M relation for the host sys- 
tems selected in both our simulations. The lines indicate 
power laws fits, whose index agrees closely with the theo- 
retical values 1/4 and 1 /3 for MONDian and Newtonian 
theory, respectively. The fitted normalisation is higher 
than the theoretical one, owning to the fact that the 
simulated halos break the hypothesis of constant density 
used in the theoretical approach. 

3.2. Normalising the relative velocities 

Special care must be taken when comparing Newtonian 
dark matter simulations to coUisionless MONDian simu- 
lations, especially when it comes to "haloes" . While we 
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Fig. 1. — Correlations of velocity dispersion (cFhoat in km/s) vs. 
effective host mass {Mj-^ast = 7.25Mhast,baryon in solar masses) 
for virialized objects of various sizes in MONDian (circles) and 
Newtonian simulations (pluses). The lines are fits with arbi- 
trary normalization for theoretically predicted scaling relations 

ai,„,f oc M^^^, , in MOND and ai,„,f oc M^^. , in 

riost,baryon iLuai, nost,baryon 

Newtonian; the best fit slopes and very close to 1/4 and 1/3 re- 
spectively. Note that the value 7.25 scales the MOND baryoninc 
mass to the Newtonian baryon plus halo mass. 

set out to use the velocity dispersion of the host system 
as the normalisation of the collision velocity in order to 
account for the mass of the host, we have just seen that 
a — M relations scales differently in MONDian than in 
Newtonian physics. What we need to do now is to move 
both simulations into the same theoretical framework 
(e.g., CDM) by appl ying the same technique outlined in 
iLlinares et al] (|2008D . i.e. we divide the MONDian (host) 
masses by the baryon fraction { ^'^^^'' = = 7.25) to 
mock dark matter haloes for direct comparison to the 
ACDM model. From these "MONDian dark matter halo 
plus baryon masses" we apply the Newtonian cr — M re- 
lation to obtain the appropriate a value to compute the 
normalized Vrei/'^- Note that we actually could reverse 
this procedure and transform the Newtonian dark matter 
haloes into the MONDian frame and use the MONDian 
a — M relation to obtain the normalized velocity. We 
have applied both methods to our data and the results 
are consistent with each other. Hence we show the results 
only for the former method to facilitate direct compari- 
son to other CDM simulations. 

4. RESULTS 

Having identified each subhalo and its corresponding 
host (cf. Section [5]), we calculate the cumulative prob- 
ability distribution of the relative velocity between host 
and its most massive subhalo Vrei ■ The result is shown in 
figure [21 We can see that the two theories have strikingly 
different behaviours. Both allow high speed satellites but 
there is a stronger tail towards high-speed in the MON- 
Dian case. The shaded lines show the high speed region, 
where the MONDian probability is about four times the 
Newtonian. 

The normalisation of the relative speed was chosen 
in order to eliminate the increased acceleration towards 
more massive host systems. The credibility of the ap- 



FlG. 2. — Accumulated probability of the relative velocity Wei 
between our host halos and their most massive subhalo normalized 
with the effective dispersion cTf^^^f = 0.01{7.25Mi-^^gf jj^j.y^„)^^^ 
km/sec (see text for explanation). 

proach detailed in Section 13.21 can now be verified by 
simply dividing our sample into different mass bins. We 
confirm that this does not lead to different results even 
though we decided to not show them in this Letter; we 
basically recover the same curves as seen in figure [2] 

Further, our results appear to be robust against slight 
changes in redshift, i.e. we neither observe a change in 
the fact that MONDian velocities are bigger nor are our 
results contaminated by the fact that we may capture 
collisions at a particular time of accidentally high ve- 
locity. The latter is confirmed by analysing the simula- 
tions at z — 0.036 leading to an indistinguishable plot. 
The same conclusion is reached when we experiment with 
other plausible normalizations or compute the distribu- 
tions of un-normalised relative speed. 

As a final test we compare our results against a Newto- 
nian model that does not contain a cosmological constant 
A, i.e. the open OCBM model of Knebe fc GibsonI ()2004f ) 
characterized by flm — 0.04 and = 0.0. We acknowl- 
edge (though not explicitly shown here) that the relative 
velocity distribution of the OCBM model is akin to the 
ACDM model presented in figure [2l we therefore ascribe 
the differences found in that plot to the effects of MOND 
rather than the (missing) cosmological constant. 

5. CONCLUSIONS 

Inspired by the observational evidence for high-speed 
encounters of galaxy clusters we studied the velocity dis- 
tribution of collisions present in two cosmological sim- 
ulations, a standard ACDM model as well as MOND. 
While there may be a problem f or ACDM t o accomodate 
such extraordinary events (e.g.. iHavashi fc White 2006) 
wc set out to quantify the probabilty for them in MOND. 
Within the limitations of our simulations, we find that 
there are substantial differences in the collision velocity 
of objects in the standard model of cosmology and its 
(possible) MONDian counterpart. Wc observe a much 
greater likelihood for high-speed collisions in MOND and 
therefore argue that this statistic can be used as a dis- 
criminator for the two competing theories. 
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We further verify numerically the velocity dispersion- 
mass relation for deep MOND gravity whose slope is dif- 
ferent to the Newtonian case {a oc M^/* for MOND in- 
stead of CT oc M^/^ for Newtonian physics). There is a 
mild scatter about these relations. 

We close with a cautionary note: the box size of 
our simulation (32ft,^^Mpc) is too small to find objects 
directly comparable to systems like the Bullet cluster. 
The coUisional velocity expected for the Bullet cluster 
{Vrei/<J ~ 2.04) is close to the limit of what we resolve 
in figure [2l While the current result is interesting, and 
the rescaled K-ei/f is likely insensitive to details of the 
simulation setup, more simulations (e.g., with a possi- 
ble neutrino component and a cosmological constant in 
MOND) are required to understand how our prediction 
depends on the cosmological model employed. 

Nevertheless, our results here may have far-reaching 
implications as well. Historically, Dark Matter and 
MOND are competing theories. Recent studies argue 
that MOND is a prescription for interactions of a cou- 
pled Dark Energy-Dark Matter field. Effectively MOND 
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is made by a non-uniform Dark Energy field. The places 
where this field condenses are identified as dark halos. 
Our results here could argue that we might differenti- 
ate between theories with interacting Dark Matter-Dark 
Energy vs. classical ACDM using data of high speed en- 
counters. It is encouraging that some subtle differences 
on how the Dark Sector self-interacts could leave signa- 
tures on "large astronomical colliders" (LAC). 
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